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Abstract 



Based upon a next-to-leading order perturbative calculation of the four-jet 
production rate in electron-positron annihilation and assuming 8 % for the 
theoretical error emerging from hadronization effects in the 0.002 < y^ut ^ 
0.004 range for the Durham clustering algorithm, we exclude the existence of 
the light gluinos at the 95 % confidence level. 
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Typeset using REVT^ 



In recent years several interesting phenomena were observed that could be naturally 
explained by the existence of light (with mass of the order of few GeV) supersymmetric 
partner of the gluon — the gluino If light gluino exists, it should be produced at 

LEP. Gluinos do not couple directly to quarks, but a gluon can decay into a gluino pair. 
Thus the numerical value of four-jet observables should be different if gluinos exist than 
predicted by pure Quantum Chromodynamics (QCD). In particular, the measurement of 
the eigenvalues of the Casimir operators of the underlying Lie group, the color charges 
is based upon measuring four-jet observables, and this should also be influenced by the 
existence of light gluinos. This fact have been utilized by all LEP experiments for setting 
exclusion limits for the existence of light gluinos However, these analyses were based 

upon leading order calculation of the four-jet observables and as a consequence the exclusion 
limits are debated [^-^l- 

Recently Dixon and Signer [§] and also ourselves have completed the calculation of 
four-jet cross sections in electron-positron annihilation at the next-to-leading order accuracy. 
In this letter we extend the second of these works to including the effect of light gluinos. 
We present the cross sections in terms of group independent kinematical functions multi- 
plying the eigenvalues of the Casimir operators of the underlying Lie group. The knowledge 
of the O(a^) corrections to the group independent kinematical functions facilitates the si- 
multaneous precision meaurement of the strong coupling and the color charge factors using 
the four-jet LEP or SLC data assuming any number of light gluino flavors. Thus in a full 
analysis the fit of the theoretical prediction to the data can reveal the number of light gluino 
flavors favoured by data. In this letter we do not perform such a full analysis. Instead, we 
point out that in the small ycnt region (?/cut < 0.004) already the four-jet rate predicted for 



the light-gluino scenario significantly differs from the recent ALEPH data [jTO|, while the 
QCD prediction is in good agreement with the experimental results. 

The main ingredients of the calculation are the four-parton next-to-leading order and 
five-parton Born level squared matrix elements. The tree level amplitudes for the pro- 
cesses e+e~ — >■ qqggg and e~^e~ — > qqQQg have been known for a long time [jll|. The 
new techniques developed by Bern, Dixon and Kosower in the calculation of one-loop mul- 
tiparton amplitudes made possible the derivation of explicit analytic expressions for 
the helicity amplitudes of the e+e^ Z^, 7* ^ 4 partons processes [|l^. Also, Campbell, 
Glover and Miller made FORTRAN programs for calculating higher order corrections to the 



e+e — i> 7* — i> 4 partons processes publicly available [0. In all of these calculations, and 



also in our work, the quark and lepton masses are set to zero. 

We use the matrix elements of refs. [|1^] for the loop corrections of the QCD subprocesses. 
As pointed out in refs. [0,^, the naive A^f = 5 — A^f = 8 shift for taking into account the 
effect of the light gluinos is not valid exactly at the next-to-leading order accuracy. Analyzing 
the color structure of the diagrams involving gluinos and using the results of refs. [l^ we 
have derived the necessary loop corrections for the process e^e~ — *• Z^, 7* — >■ qqgg and 
the modifications induced by the presence of the gluinos on the four-parton one-loop QCD 
amplitudes |]T^|. We have also derived the e+e~ qqggg tree-level amplitudes. 



It is well known that the next-to-leading order correction is a sum of two integrals — 
the real and virtual corrections — that are separately divergent (in the infrared) in d = 4 
dimensions. For infrared safe observables, for instance the four-jet rates used in this work, 
their sum is finite. In order to obtain this finite correction, we use a slightly modified version 
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of the dipole method of Catani and Seymour [|T6| that exposes the cancellation of the infrared 
singularities directly at the integrand level. The formal result of this cancellation is that the 
next-to-leading order correction is a sum of two finite integrals, 
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where the first term is an integral over the available five-parton phase space (as defined by 
the jet observable) and the second one is an integral over the four-parton phase space. 

Once the phase space integrations in eq. (|l|) are carried out, the next-to- leading order 
differential cross section for the four-jet observable O4 takes the general form 
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In this equation ctq denotes the Born cross section for the process e+e^ qq, s is the total 
cm. energy squared, /i is the renormahzation scale, while and Co, are scale independent 
functions, i^o, is the Born approximation and Co, is the radiative correction. We use the 
two-loop expression for the running coupling. 
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Pi = - ^CpTnNi - -CaTrN, - -C\N,^ , (6) 

for number of light gluino flavors [0, and with the normalization Tr = 1/2 in 
Tr(T"T''"'') = Tpid""^ . The numerical values presented in this letter were obtained at the 

peak with Mz = 91.187 GeV, Tz = 2.49 GeV, sin^^ = 0.23, a,(Mz) = 0.118 and 
A^f = 5 light quark flavors. 

In the case of jet rates it is customary to normalize the cross section to the 0{as) total 
cross section, cxtot = ctq (1 + Q^s/ti")- Then the four-jet rate can be written as 
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where the scale dependence in as has been suppressed. The Born approximation and the 
higher order correction are hnear and quadratic forms of ratios of the color charges ||15|| : 



B^ = Bo + B^x + ByV , (8) 

and 

C^ = Co+ C^x + Cyy + C,z (9) 
+ Cxx X + Cxy xy -\- Cyy y . 

The X and y parameters are ratios of the quadratic Casimirs, x = Ca/Cf and y = Tr/Cf, 
while z is related to the square of a cubic Casimir, 

Na 

C^= J2 Tr(T"T''T"^")Tr(r^'=T''T") , (10) 

a,b,c=l 

via z = TTr^r- We plot the numerical values for the correction functions Ci for the case 



of Durham clustering algorithm |]T8[ in Fig. 1. Note that the Cz function is completely 
negligible, therefore it is not plotted. In order to exhibit the behaviour in the small T/cut 
region, the correction functions are plotted against L^, where L = ln(l/?/cut)- 

Using these group independent kinematical functions one can fit the observed data for 
the best values of the color factor ratios x and y in both theories separately, and test 
how the obtained results are compatible with the SU(3) values x = 9/4 and y = 3/8. 
In principle, any kind of four-jet observable can be used for this purpose. In practice, 



various normalized angular distributions were used These observables have the virtue 
that their renormalization scale dependence is very small already at tree level due to the 
normalization, therefore they can be calculated reliably at leading order. Indeed, this has 
been shown to be the case in the next-to-leading order calculation in ref. pO| and we confirm 
that result. We should like to emphasize however, that the small scale dependence of the 
tree-level perturbative result does not necessary mean that the results for the exclusion of 
light gluinos from color charge measurement are insensitive to the higher order corrections. 
Fig. 1 shows that the C^, C^x and Cxy kinematical functions change substantially, while 
the Cy and Cyy functions remain unchanged with the inclusion of the additional degrees of 
freedom. This may modify the result of the fit for the color charges. 

We do not perform the complete color charge analysis here, but note that this approach 
is slightly different from the effective number of flavors approach, where the same best fit 
for both theories is made and this unique result is compared with x = 9/4 and y^Q = 3/5 
with the effective ?/efr value obtained from the A^f = 5 — >■ 8 rule. This latter method is made 
possible by the simplicity of the tree-level matrix elements, but strictly speaking is not valid 
in a next-to-leading order analysis. 

In this letter we simply take values in SU(3) for the color factor ratios and using our 
kinematical functions we predict the numerical values for the Born and correction functions 
BAiycnt) and C^iycut) in both theories. These functions are then fitted according to the 



formula 22 



Fivcut) = E > (11) 

n=0 
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where F = B4, or C4. The coefficients fc„ are tabulated in Table I. The fit-range is 0.002-0.1. 
The QCD results agree with the rates obtained by Dixon and Signer [^. 

Using the B4 and C4 functions, we construct Ri^i/cut) according to formula (0). We see 
that for small values of resolution parameter, ?/cut < 0.004, the predictions for the QCD 
and QCD + light gluino case start to differ. For small values of j/cut however, the fixed 
order perturbative prediction is not reliable, because the expansion parameter as In^ ?/cut 
logarithmically enhances the higher order corrections. The usual remedy is the all order 
resummation of the leading and next-to-leading logarithmic contributions and matching the 
resummed and fixed order results |^T|. We use R- matching as described in refs. [^1],^] . 

Fig. 2 shows the effect of resummation and matching in QCD (A^ig = 0). We see that 
the resummed and matched prediction agrees with the ALEPH data corrected to hadron 
level []T0[ within hadronization uncertainty. Encouraged by this, in the same figure we plot 
the fixed order and the resummed and matched SU(3) predictions for the QCD + light 
gluino scenario (A^ig = 1). The bands indicate the theoretical uncertainty of the fixed order 
prediction due to the variation of the renormalization scale = /i/i/s between 1/2 and 
2. The data are well described by the QCD prediction in the whole Ucut range, but they 
fall outside the QCD + light gluino band. For values of |/cut < 0.004 the theoretical bands 
bifurcate making the discrepancy between the two theories significant. 

In order to quantify the exclusion level of the very light gluino, we plot the scale depen- 
dence of the next-to-leading order prediction at a fix value of ?/cut = 0.002 in Fig. 3. Our 
choice for this value of the resolution parameter comes from the observation that in Fig. 2 
the intercept of the iVig = theoretical curves for the O(as^) and 0(as^)+NLLA prediction 
at = 1 is around 0.002, showing that in this region the fixed order prediction itself can 



be trusted. Also around ?/cut = 0.002 the hadronization effects are still small p3| , p4 . 

From Fig. 3 we see that the QCD prediction agrees well with the measured data, while 
the QCD + light gluino prediction falls significantly below. Taking the renormalization 
scale dependence as the theoretical error, the light gluino is excluded at more than 4(j, or 
99.99 % confidence level. The theoretical prediction for the two theories differ at the ~ 85 % 
confidence level. The j/cut = 0.002 value at LEP corresponds to a fc^ ~ 5 GeV, where the 
theoretical error emerging form hadronization does not exceed 8% [Q. Therefore, with 
hadronization uncertainty included the QCD + light gluino prediction is still excluded at 
the 2a, or 95 % confidence level. According to Fig. 2, the effect of resummation of leading 
and next-to-leading ln(l/?/cut) logarithms is to make the exclusion even more significant. 

Finally, we would like to remark that the inclusion of the non-zero gluino mass decreases 
the tree-level cross section [0]. The K factor for the QCD -|- light gluino theory is small, 
K ^ 1.25, suggesting that this mass effect is bound to be robust when including the higher 
orders. 

In this letter we presented the next-to-leading order corrections to the gauge indepen- 
dent kinamatical functions of the four-jet rate cross section in QCD and in QCD with the 
extension of light gluinos for the case of the Durham clustering algorithm. The knowledge 
of these corrections facilitates the simultaneous precision meaurement of the strong coupling 
and the color charge factors using the four-jet LEP or SLC data assuming any number of 
light gluino flavors. Assuming standard SU(3) values for the color factors, we presented the 
four-jet rates as predicted by next-to-leading order perturbation theory matched with all 
order resummed results in the next-to-leading logarithmic approximation. These theoreti- 
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cal predictions were compared to recent data obtained by the ALEPH collaboration. The 
four-jet rates in the QCD and QCD + light ghiino scenario are significantly different for 
small values of the resolution parameter i/cut- The data favour the gluinoless world, exclud- 
ing the existence of light gluinos at the 95 % level. This result was argued to be robust 
against hadronization and gluino mass effects. Our results also show that a conventional 
color charge measurement based on angle distributions should give the most significant light 
gluino exclusion limit if the jets are defined at small i/cutj (^cut < 0.004). 

We thank L. Dixon for providing Mathematica files of the one-loop helicity amplitudes, 
furthermore S. Catani and G. Dissertori for helpful correspondence. This research was sup- 
ported in part by the EEC Programme " Human Capital and Mobility" , Network " Physics at 
High Energy Colliders", contract PECO ERBCIPDCT 94 0613 as well as by the Hungarian 
Scientific Research Fund grant OTKA T-016613. 
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FIG. 1. Group independent kinematical functions Cj for QCD and QCD + light gluino theory. 



TABLES 



TABLE I. CoefHcients of the analytic fits for the Born level and next-to-leading order scale 
independent functions B4 and C4. First two rows for QCD, second two rows for QCD -|- light 
gluino theory. 
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FIG. 2. Full next-to-leading order and resummed-matched four-jet rate in QCD -|- Nig = 0, 1 
compared to ALEPH data. 
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FIG. 3. Renormalization scale dependence of the next-toleading order prediction for the four-jet 
rate at ycut = 0.002. The grey band indicates the value measured by ALEPH, corrected to hadron 
level. 
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